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Abstract — In 2012, during the Entry, Descent, and Landing 
(EDL) of the Mars Science Laboratory (MSL) entry vehi- 
cle, a 21.5 m Viking-heritage, Disk-Gap-Band, supersonic 
parachute will be deployed at approximately Mach 2. The 
baseline algorithm for commanding this parachute deploy- 
ment is a navigated planet-relative velocity trigger. This pa- 
per compares the performance of an alternative range-to-go 
trigger (sometimes referred to as ’’Smart Chute”), which can 
significantly reduce the landing footprint size. Numerical 
Monte Carlo results, predicted by the POST2 MSL POST 
End-to-End EDL simulation, are corroborated and explained 
by applying propagation of uncertainty methods to develop an 
analytic estimate for the standard deviation of Mach number. 
A negative correlation is shown to exist between the standard 
deviations of wind velocity and the planet-relative velocity at 
parachute deploy, which mitigates the Mach number rise in 
the case of the range trigger. 
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Nomenclature 

fi statistical mean 

p atmospheric density 

p xy correlation coefficient 

<j standard deviation 

a 1 2 statistical variance 

c s speed of sound 

M Mach number 

q dynamic pressure 

Va atmospheric-relative velocity 

V R planet-relative velocity 



Figure 1. Artist concept of Mars Science Laboratory rover, 
named Curiosity , exploring Mars. 


Vw 

wind velocity 

DGB 

Disk-Gap-Band Parachute 

EDL 

Entry, Descent, and Landing 

IMU 

Inertial Measuring Unit 

MER 

Mars Exploration Rover 

MOLA 

Mars Orbiting Laser Altimeter 

MSL 

Mars Science Laboratory 

POST2 

Program to Optimize Simulated Trajectories II 

TPS 

Thermal Protection System 


1. Introduction and Background 

Mars Science Laboratory 

Scheduled to launch in the 2011 Type-I Mars opportunity, 
Mars Science Laboratory (MSL) will be the next robotic mis- 
sion to explore the surface of the red planet under NASA’s 
Mars Exploration Program. Building on the success of the 
twin Mars Exploration Rover ( MER ) rovers, Spirit and Op- 
portunity , which landed in 2004, the MSL rover, named Cu- 
riosity. , will assess the planet’s present and past habitability. 
Armed with a suite of state-of-the-art scientific instruments, 
Curiosity , shown in Figure 1 , will collect dozens of soil and 
rock samples and analyze them in an on-board laboratory for 
the presence of biomarkers, the chemical building blocks of 
life. 


1 U.S. Government work not protected by U.S. copyright. 

2 IEEEAC Paper #1 142, Final Version 1, Updated January 5, 2011. 
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Designing an Entry, Descent, and Landing (EDL) system to 




Figure 2. Mars Science Laboratory DGB parachute under- 
going full-scale wind-tunnel testing. 


deliver such a large and capable rover safely to a scientifi- 
cally compelling site, which is rich in minerals likely to trap 
and preserve biomarkers, presents a myriad of engineering 
challenges. Not only is the payload mass significantly larger 
than all previous Mars missions, the delivery accuracy and 
terrain requirements are also more stringent. In August of 
2012, MSL will enter the Martian atmosphere with the largest 
aeroshell ever flown to Mars, fly the first guided lifting entry 
at Mars, generate a higher hypersonic lift-to-drag ratio than 
any previous Mars mission, and decelerate behind the largest 
supersonic parachute ever deployed at Mars. The MSL EDL 
system will also, for the first time ever, softly land Curiosity 
directly on her wheels, ready to explore the planet’s surface. 

Parachute Decelerators for Mars 

Since the first Viking landing in 1976, the super-sonic deploy- 
ment of a parachute has been a critical event in all Mars EDL 
systems. This is because at Mars, due to the planet’s thin 
atmosphere, only entry systems with ballistic coefficients be- 
low about 50 kg/m 2 have the ability to deliver payloads to 
subsonic terminal velocities [1]. For MSL , a blunt aeroshell, 
with a ballistic coefficient of approximately 140 kg/m 2 , is 
first used to slow the vehicle from hypersonic entry veloci- 
ties as high as 6 km/ s down to low super-sonic speeds, near 
400 m/s. At that point, the 21.5 m diameter Disk-Gap-Band 
Parachute Decelerator System (PDS), shown in Figure 2, is 
then used to reduce the ballistic coefficient to approximately 
15 kg/m 2 . The parachute continues slowing the vehicle be- 
low Mach 1 to a sub-sonic terminal velocity of approximately 
100 m/s. 

Because of the importance of the parachute deployment 
event, parachute failure is a key risk considered in Mars EDL 
system design. Higher Mach numbers and dynamic pres- 
sure during parachute inflation put the parachute at a higher 
risk of failure due to three factors: (1) higher dynamic pres- 
sures result in higher structural loads on the parachute; (2) 


higher Mach numbers result in increased aerothermal heating 
of parachute structure, which can reduce material strength; 
and (3) at Mach numbers above Mach 1.5, DGB parachutes 
exhibit an instability, known as areal oscillations, which re- 
sult in multiple partial collapses and violent re-inflations. 
The chief concern with high Mach number deployments, for 
parachute deployments in regions where the heating is not a 
driving factor, is therefore, the increased exposure to areal 
oscillations. 

The Viking parachute system was qualified to deploy between 
Mach 1.4 and 2.1, and a dynamic pressure between 250 and 
700 Pa [1]. However, Mach 2.1 is not a hard limit for suc- 
cessfully operating DBG parachutes at Mars and there is very 
little flight test data above Mach 2.1 with which to quantify 
the amount of increased EDL system risk. Figure 3 shows 
the relevant flight tests and flight experience in the region of 
the planned MSL parachute deploy. While parachute experts 
agree that higher Mach numbers result in a higher probabil- 
ity of failure, they have different opinions on where the limit 
should be placed. For example, Gillis [5] has proposed an up- 
per bound of Mach 2 for parachute aerodynamic decelerators 
at Mars. However, Cruz [3] places the upper Mach number 
range somewhere between two and three. 

This presents a challenge for EDL system designers, who 
must then weigh the system performance gains and risks as- 
sociated with deploying the parachute earlier, at both higher 
altitudes and Mach numbers, against a very real, but not well 
quantified, probability of parachute failure. It is clear that 
deploying a DGB at Mach 2.5 or 3.0 represents a significant 
increase in risk over an inflation at Mach 2.0. However, it is 
not clear how much additional risk is encumbered by deploy- 
ing the parachute at 2.25 instead of 2.05. This is especially 
true for Mars EDL in light of the extremely large uncertainties 
in the flight environment, especially atmospheric density and 
winds, that result in very large uncertainties in Mach number. 

Parachute Deployment Algorithms 

Previous missions have utilized various methods for trigger- 
ing parachute deployment. Though the parachute qualifica- 
tion has been stated in terms of Mach number and dynamic 
pressure, no previous mission has had the ability to directly 
measure either of these quantities. Therefore, all missions 
have had to rely on proxy measurements of other states in or- 
der to infer whether or not conditions were safe for deploying 
the parachute. Viking used a radar altimeter measurement to 
trigger this critical event [6]. Mars Pathfinder [8] and MER 
[4] both used triggers based on sensed-acceleration measure- 
ments, provided by the on-board Inertial Measurement Unit 
(IMU), to proxy dynamic pressure, though their algorithms 
differed. Mars Phoenix Lander , on the other hand, used a 
navigated velocity trigger to proxy Mach number. 

When originally proposed, MSL (known then as Mars Smart 
Lander) featured a range trigger as part of the Apollo - heritage 
entry guidance system, often referred to as the Smart Chute 
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Table 1. MSL Candidate Landing Sites 
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Figure 3. Relevant test and flight experience of supersonic 
Disk-Gap-Band parachutes in the region of MSL parachute 
deployment (shaded region). 


algorithm. To minimize the parachute deploy footprint, the 
parachute was to be deployed on a command from the en- 
try guidance when the estimated range-to-go to the target 
was minimized [7]. As an additional safety measure, the 
Smart Chute algorithm also included navigated velocity lim- 
its. These limits were used in order to protect the parachute 
against either excessively high dynamic pressures or exces- 
sively low deployment altitudes [2]. Above the high veloc- 
ity set-point, parachute deploy was inhibited. Below the low 
velocity limit, parachute deploy was triggered, regardless of 
range-to-go. 

Eventually, the Smart Chute range trigger was dropped from 
the MSL baseline in-favor of a velocity trigger. The rationale 
for this decision was to maximize the altitude performance 
of the system, which was being strained at that time by rapid 
mass growth of the rover and a very challenging altitude re- 
quirement for demonstrating the capability to land as high as 
+2.0 km above the MOLA reference areoid. It was argued 
at the time, that due the monotonically decreasing altitude 
and velocity just prior to parachute deploy, the upper velocity 
limit of the Smart Chute represented the earliest, and there- 
fore highest, deployment condition that was considered safe. 
Replacing the Smart Chute trigger with a pure velocity trig- 
ger, operated at the same set-point as the upper velocity limit, 
would maximize the parachute deploy altitude, while main- 
taining the same level of risk to the parachute. 

Study Motivation 


Site 

Name 


Mawrth Valis 
Gale Crater 
Eberswalde Crater 
Holden Crater 


Lat. 

(deg) 


24.01°N 
4A9°S 
23.86 °S 
26.37 °S 


Lon. 

(deg) 


341.03°E 
137A2°E 
326.73 °E 
325.10 °E 


Elevation 

(km) 


-2.25 

-4.45 

-1.45 

-1.94 


900 


pose of proposing and selecting possible landing sites. While 
many sites were initially proposed at the first of these work- 
shops, the outcome of the 4th Landing Site Workshop in 2008 
was a list of four candidate sites, listed in Table 1. Of the 
four final sites, Eberswalde Crater has the highest elevation 
at -1.45 km MOLA , which is significantly below the altitude 
capability of the system (estimated to be somewhere around 
0 km MOLA). These lower site altitudes have improved EDL 
timeline margins significantly compared to the time when the 
parachute deploy trigger was changed. In light of this reduced 
premium on altitude, this study sought to re-evaluate the mer- 
its of a range trigger relative to the baseline velocity trigger. 

2. Monte Carlo Results 

A side-by-side comparison of a range trigger and velocity 
trigger was conducted for MSL. The purpose of this com- 
parison was to evaluate the relative performance of the range 
trigger relative to the baseline velocity trigger. A single hy- 
brid velocity-range trigger was developed that is capable of 
emulating either velocity or range triggers by appropriate 
choice of parameters. This trigger works by specifying a lin- 
ear switching curve in velocity-range space. The algorithm 
triggers parachute deploy when the velocity drops below the 
switching curve for the given range to target. A horizontal 
switching curve, therefore, produces a pure velocity trigger, 
while a vertical switching curve, on the other hand, produces 
a pure range trigger. 

For each trigger a 6-DoF Monte Carlo analysis was per- 
formed using the 08-GAL-06 MSL POST2 end-to-end EDL 
performance simulation. The two triggers were each indepen- 
dently tuned to produce the same nominal parachute deploy at 
Mach 2.0, as was the standard project procedure for running 
Monte Carlos. It was expected that the results would show 
a smaller parachute deploy footprint for the range trigger at 
the expense of reduced altitude performance and increased 
deploy Mach number. 


As stated, the switch from a range trigger to a velocity trigger 
had been argued for the maximization of parachute deploy 
altitude. Though the project would eventually receive some 
relief from the +2.0 km altitude requirement, a premium on 
altitude performance existed for quite some time. Concur- 
rently, however, the project had initiated a series of Landing 
Site Workshops, open to the scientific community, for the pur- 


Figure 4 shows the expected reduction in the 99.5%-tile foot- 
print ellipse. In this case, the ellipse was reduced from 16.7 
by 7.5 km for the velocity trigger to 7.7 by 4.1 km for the 
range trigger, a 75% reduction in area. However, the expected 
altitude loss and Mach increase were not observed. 

Table 2 contains a summary of Monte Carlo results for the 
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Velocity Trigger 


Table 2. Monte Carlo Results 



^ i i i i i i i i i i 

22 20 18 16 14 12 10 8 6 4 2 0 

Down-range to Target (km) 

Range Trigger 



^ i i i i i i i i i i 

22 20 18 16 14 12 10 8 6 4 2 0 


Down-range to Target (km) 

Figure 4. Comparison of Parachute Deploy Footprints. 


velocity and range triggers. The mean values in Table 2 are 
similar, by design, because both triggers were tuned for the 
same nominal parachute deploy Mach number. However, 
even though the standard deviation of planet-relative veloc- 
ity is nearly 5.5 times larger for the range trigger, the stan- 
dard deviation of atmospheric-relative velocity is only 24% 
higher, and the standard deviation of Mach number increased 
by only 0.006. This results in an increase of only 0.033 in the 
99th percentile of Mach number. Figure 5 compares the his- 
tograms for parachute deploy Mach number. Strictly speak- 
ing, even this slight increase in Mach number represents some 
higher risk to the parachute system. Practically, however, 
the amount of additional risk is negligible and impossible to 
quantify given the limited test data available. Therefore, the 
nominal deploy Mach number for the range trigger was not 
lowered to bring the 99% -tile Mach in agreement with the 
velocity trigger results. 

Finally, because the mean deploy altitudes were identical 
(within 10 m) and because the standard deviation of altitude 
was slightly smaller for the range trigger, the range trigger 
showed a 218 m increase in 1%-tile altitude over the velocity 
trigger. This result was contrary to the expected altitude loss, 
but related to the negligible increase in Mach number. An 


Parameter 

Units 

Velocity 

Trigger 

Range 

Trigger 

Mean Value 

Planet-rel. velocity 

m/s 

463.8 

464.2 

Atm. -relative velocity 

m/s 

447.2 

447.3 

Mach number 

- 

2.002 

2.001 

Dynamic pressure 

Pa 

546.7 

546.7 

Parachute opening load 

1000 lbs. 

47.64 

47.67 

Altitude wrt. MOLA 

km 

10.084 

10.073 

Standard Deviation 

Planet-rel. velocity 

m/s 

3.0 

16.2 

Atm. -relative velocity 

m/s 

13.0 

16.1 

Mach number 

- 

0.096 

0.102 

Dynamic pressure 

Pa 

37.2 

39.5 

Parachute opening load 

1000 lbs. 

3.57 

3.43 

Altitude wrt. MOLA 

km 

1.129 

1.077 

99 th Percentile 

Planet-rel. velocity 

m/s 

470.6 

499.1 

Atm. -relative velocity 

m/s 

475.1 

483.6 

Mach number 

- 

2.221 

2.254 

Dynamic pressure 

Pa 

639.9 

647.1 

Parachute opening load 

1000 lbs. 

56.92 

56.28 

1 st Percentile 

Altitude wrt. MOLA 

km 

7.343 

7.561 


appreciable increase in 99% -tile Mach number would have 
forced a reduction in mean Mach number, resulting in a mean 
altitude loss of approximately 85 m per 0.01 Mach number 
decrease. 

These Monte Carlo results suggested a very favorable trade to 
adopt the range trigger, but more analysis was needed to ex- 
plain why they were contrary to the expectations of the EDL 
team. However, because Monte Carlo is a numerical analy- 
sis, it is very difficult, if not impossible, to explain why these 
results differed so dramatically from expectation. Therefore, 
an analytical approach is developed in the following section 
to attempt to explain the observed numerical results. 

3. Theoretical Development 

Analytical Approach 

We desire an analytical expression which relates random vari- 
ables at the instant of parachute deploy to the resulting distri- 
bution in Mach number. To simplify the analysis, we will 
assume that all random variables are normally distributed. 
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Velocity Trigger 



Mach Number 



Mach Number 

Figure 5. Comparison of Mach Number at Parachute De- 
ploy. 

Though this assumption is not correct in general, the partic- 
ular random variables of interest from the numerical Monte 
Carlo simulation are approximated very well by Gaussian 
distributions. Furthermore, we will focus on estimating the 
standard deviation, rather than the mean, since the later is a 
tunable parameter. 

We wish to write, therefore, the standard deviation of Mach 
number, cjm, as a function of the standard deviations of rel- 
ative velocity ( Vr ), wind velocity (Vw), and speed of sound 
(c s ), ay R , °v W ' an d &c s respectively. We see from the scalar 
definition of atmospheric-relative velocity (Eqn. 1) and Mach 
number (Eqn. 2) that these flight parameters are the con- 
stituents of Mach number. Here V a is defined as a scalar sum, 
rather than a strict vector summation, so that Vw is more cor- 
rectly the projection of the wind vector on the planet-relative 
velocity vector, with the convention of a positive value indi- 
cating a head-wind. 


V A = V r + Vw 

(1) 

§ 

II 

■1^ 

(2) 




Propagation of Uncertainty 

We apply propagation of uncertainty to the above definitions 
to get expressions for the standard deviations of atmospheric- 
relative velocity and Mach number. We assume that distri- 
butions of the flight conditions, Vr, Vw, and c s , are either 


known or specified and that they are normally distributed, 
such that V R = N (fi VR ,cry R )> V w = N (ji Vw , ct ^ w ) , and 
c s = N (fi Cs ,a^). 

Note that two of these random variables, wind velocity and 
speed of sound, are environmental parameters. Their distri- 
butions are outputs of the simulation’s atmosphere model in 
the region of the parachute deploy point and are independent 
of the trigger algorithm. For MSL , the uncertainty in these pa- 
rameters is determined statistically by sampling multiple sols 
and local solar times predicted by meso- scale computer simu- 
lations of the Martian atmosphere. The third random variable, 
planet-relative velocity, is a result of integrating the equa- 
tions of motion in time from the initial state to the instant 
of parachute deploy. Since different algorithms will trigger 
parachute deploy at different times, the uncertainty in Vr is a 
characteristic of the algorithm employed. 

Equation 1 may be recognized as a linear combination of nor- 
mal distributions. Therefore, the result is known to be also 
normally distributed, Va = N (py A , a y A )’ with a standard 
deviation ay A given by Equation 3. Note that while the dis- 
tributions Vr and Vw are assumed to be normal distributions, 
they are not necessarily independent. The angle, 0, in the 
third term of the radical is a correlation angle, which is re- 
lated to the correlation coefficient, p xy , between the two ran- 
dom variables. A correlation angle of 6 = 0 is equivalent to a 
correlation coefficient of p xy = 0, which implies that there is 
no linear correlation between the variables. Likewise, a corre- 
lation angle, 6 = | , is equivalent to a correlation coefficient 
of Pxy = 1, which implies a direct one-to-one relationship 
between the variables. 

°v A = f <4 r + Vw + VWv w sin 0 (3) 

Similarly, by taking the natural logarithm of Equation 2, 
Equation 4 may also be recognized as a linear combination 
of normal distributions, if one assumes that the logarithm of 
a normal distribution is also Gaussian. This is approximately 
true if the random variable has a non-zero mean and a stan- 
dard deviation that is small compared to that mean. Equation 
5 relates the standard deviation of the natural logarithm of a 
random variable to the random variable’s standard deviation. 
In practice, the error in this approximation is typically less 
than 1% for random variables that have a mean to standard 
deviation ratio (px/&x) greater than about 12. 

In M = In Va — In c s (4) 


CTlnW w W\ ax (5) 

Using this approximation results in Equation 6 for the stan- 
dard deviation of In M. Again, the random variables are all 
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assumed to be normally distributed, but not necessarily un- 
correlated. The angle, a , in the third term of the radical in 
Eqn. 6 is the correlation angle between In Va and In c s . 


O'ln M = y crpn y A + of n Cs + 2 y A a ln Cs sin a (6) 

Finally, substituting Eqn. 3 into Eqn. 6 and using the ap- 
proximation, Eqn 5, results in the desired expression for <tm- 
This expression, Equation 7, predicts the standard deviation 
of Mach number as a function of the three mean values (pv R , 
Pv w , and /i Cs ), three standard deviations (cry R , cry w , and 
<r Cg ), and the two correlation angles (0 and a). 



+B sin a 

a = 2 <7Vr<7Vw 

(Mv r + MVw ) 

B 2 (^k + + 2 fv^Vw sin d ) * g Ca 

Vc s 

The first term in the radical of Eqn. 7 is the planet-relative ve- 
locity term. The magnitude of the variance of Vr is not only 
dependent on the integrated effect of the dispersions leading 
up to parachute deploy, but also on the particular triggering 
algorithm, since these will affect the time at which deploy- 
ment occurs. The second and third terms are the wind ve- 
locity term and speed-of- sound term, respectively, which are 
independent of the triggering logic. The fourth term is the 
wind correlation, which will be discussed in more detail in 
the following section. This is the statistical correlation be- 
tween the wind velocity and the planet-relative velocity at the 
moment of parachute deploy. The last term is the speed-of- 
sound correlation. 

4. Discussion 

Comparison of Results 

The analytic expressions derived in the previous sections 
were validated by comparison with the Monte Carlo numeri- 
cal results. Three Monte Carlo results are compared: a veloc- 
ity trigger, a hybrid trigger, and a range trigger. For the hybrid 
trigger case, the switching curve was set with a finite, positive 
slope to produce intermediate results that were neither a pure 
velocity trigger nor a pure range trigger, but roughly half- 
way between the two. Figure 6 plots the resulting dispersions 
in the velocity-range plane and shows the primary effect of 
varying the trigger, which is to slice the corridor of incom- 
ing trajectories at a different angle changing the projections 



Figure 6. Parachute Deploy Conditions, Velocity vs. Range. 

of the dispersions in the altitude and velocity axes. 

Statistics from the Monte Carlos, presented in Table 3, were 
collected on flight parameters Vr, Vw, and c s at parachute 
deploy as well as the two correlation angles, 6 and a. These 
correlation angles are found from the arcsine of the correla- 
tion coefficients of the Monte Carlo data. These values were 
used to estimate a m (Equation 7). As shown in Table 4, these 
analytic estimates compared very well with the numerical re- 
sults with a maximum error of 0.3% of the actual result. 

Table 3 shows that <jy R progressively increases from 3.0 m/s 
for the velocity trigger to 16.2 m/ s for the range trigger. The 
wind dispersions, ay w are consistent, as expected, since it is 
an environmental parameter. The slight decrease in a Cs is due 
to altitude differences. Note that the wind correlation angle, 
6, is negative and that its value is the same for the hybrid and 
range triggers and only a few degrees different for the velocity 
trigger. This correlation is discussed in the following section. 

Table 3. Parachute Deploy Flight Conditions 


Parameter 

Units 

Velocity 

Trigger 

Hybrid 

Trigger 

Range 

Trigger 

Mean Value 

Planet-rel. vel. 

m/s 

463.8 

463.7 

464.2 

Wind velocity 

m/s 

-16.6 

-16.8 

-16.9 

Speed-of-sound 

m/s 

223.7 

223.8 

223.8 

Standard Deviation 

Planet-rel. vel. 

m/s 

3.0 

8.5 

16.2 

Wind velocity 

m/s 

13.6 

13.7 

13.7 

Speed-of-sound 

m/s 

8.0 

7.9 

7.8 

Correlation 

o (Vr, V w ) 

deg 

-18.0 

-25.6 

-25.4 

a (V A , c s ) 

deg 

4.6 

3.2 

1.7 
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Table 4. Comparison of Analytic and Numeric Results for 
Standard Deviation of Mach Number 


& M 

Analytic 

Numerical 


Estimate 

Result 

Velocity Trigger 

0.0957 

0.0959 

Hybrid Trigger 

0.0931 

0.0933 

Range Trigger 

0.1019 

0.1022 


Wind Correlation 

Because the velocity trigger has, by design, a negligible stan- 
dard deviation of planet-relative velocity (the triggering pa- 
rameter), the variance in atmospheric-relative velocity is due 
solely to the variance in the wind at parachute deploy. In con- 
trast, the range trigger has a much larger variance in planet- 
relative velocity, which is traded, again by design, for a small 
variance in range. Consequently, the range trigger, having the 
same wind dispersions as the velocity trigger, should be ex- 
pected to have a much larger variance in atmospheric-relative 
velocity, and therefore, Mach number. However, this result is 
not seen in the Monte Carlo analysis. 

The reason for this is the negative correlation angle ( 6 ) be- 
tween the planet-relative velocity and the wind velocity at 
parachute deploy, shown in Table 3. Physically, this corre- 
lation exists because the drag, experienced during entry, is a 
function of the atmospheric-relative velocity, and therefore, 
wind. If the entry vehicle is experiencing a head wind (posi- 
tive wind speed by our convention), the atmospheric-relative 
velocity and drag will be higher than a no- wind case. There- 
fore, the head-wind case will arrive at the parachute deploy 
point at a lower planet-relative velocity. Conversely, if the 
entry vehicle experiences a tail wind (negative wind), the 
atmospheric-relative velocity and drag will be lower, result- 
ing in a higher planet-relative velocity at parachute deploy. 

Graphically, the effect of this wind correlation on ay A can be 
seen in the example plotted in Figure 7. Here the blue vector 
(OA) has been drawn along the horizontal axis with a magni- 
tude equal to the standard deviation of the wind speed, ay w 
- 15 m/s, in this example. The standard deviation of planet- 
relative velocity (OB) is drawn in green, ay R = 20 m/s , 6 
= -22 deg. The angle between vectors OA and OB is deter- 
mined by the correlation between the random variables Vr 
and Vw- The red vector (OC) is the standard deviation of 
atmospheric-relative velocity, ay A =20 m/s, resulting from 
the vector addition of OA and OB. Because of the negative 
correlation, ay A is significantly smaller than the 25 m/s that 
would be expected for uncorrelated random variables. 

For all hybrid triggers, the magnitude of ay R is smaller, but 
the correlation angle, 6 is the same. This is because the cor- 
relation angle is a property of the simulation, independent of 
the trigger, resulting from the integration of the winds prior to 
parachute deploy. Therefore, the chord AC in Figure 7 repre- 


0 = 0 * 



Figure 7. Polar representation of Equation 3 plotted for a 
range trigger. Standard deviation of wind speed (OA) and 
standard deviation of relative velocity (OB) add vectorially. 
The result is standard deviation of atmospheric velocity (OC). 


sents the locus of all possible hybrid trigger points, C, as the 
magnitude of vector OB is reduced from 20 m/s (the range 
trigger limit) to zero (the velocity trigger limit). In this ex- 
ample, for hybrid triggers with ay R less than 11.2 m/s, the 
resulting ay A is smaller than the velocity trigger result. Sim- 
ilarly, the minimum oy A is achieved when ay R is approxi- 
mately 5.6 m/ s. 

Risk Assessment 

In the previous section the negative wind correlation was 
shown to mitigate the Mach number increase for range trig- 
gers. This explains and validates the numerical Monte Carlo 
results. However, the presence of this correlation creates ad- 
ditional considerations in determining the risk posture of the 
EDL system. While both velocity and range triggers are af- 
fected by dispersions in the wind model, the manner in which 
they are affected is different. The velocity trigger, by means 
of a near-zero ay R , is affected only by the instantaneous wind 
at the moment of parachute deploy. Conversely, the range 
trigger, by means of the wind correlation term, is affected by 
the integrated history of the winds throughout the entire entry. 

This levies stricter requirements on the wind modeling and 
places an additional lien on simulation results. What if the 
correlation angle is not realistic or not even present at all? In 
a worst case, an uncorrelated wind (6 = 0 deg) would increase 
the estimate of gm in Table 4 from 0.1019 to 0.1178, rais- 
ing the estimated 99% -tile Mach number by approximately 
0.05. Using the previous thumb-rule for altitude, this could be 
countered by an approximate 435 m decrease in nominal alti- 
tude. However, this retuning of the nominal must be weighed 
against the risk of parachute failure if the Mach increase in 
simply accepted. Returning to Figure 3, there is little evi- 
dence that a 0.05 increase in Mach number would increase 
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risk appreciably. 

5. Summary 

A side-by-side comparison of parachute deployment triggers 
was conducted for MSL. Two triggers, the baseline velocity 
trigger and an alternate range trigger, were tuned to produce 
the same nominal parachute deploy Mach number. For each 
trigger a 6-DoF Monte Carlo analysis was performed using 
the MSL POST2 end-to-end EDL performance simulation. 
The results of this study showed that a range trigger has the 
potential to significantly reduce footprint size with negligible 
Mach increase and no altitude loss. 

An expression was developed, Equation 7, which predicts the 
standard deviation of Mach number, gm as a function of the 
three mean values (p,y R9 pty W9 and fi Cs ), three standard de- 
viations (Gy R , Gy w , and g C s ), and two correlation angles (6 
and a). This expression was shown to have very good agree- 
ment with the numerical Monte Carlo results and was used 
to explain why Mach number failed to increase for the range 
trigger. 

It was shown that a negative correlation coefficient exists 
between the planet-relative velocity and the wind velocity 
at parachute deploy. Physically, this correlation exists be- 
cause the drag, experienced during entry, is a function of the 
atmospheric-relative velocity, which includes the wind. In 
the case of a range trigger, this wind correlation mitigates 
the rise in Mach number uncertainty due to the variance of 
planet-relative velocity inherent to the trigger. 

Finally, Equation 7 was again used to estimate the additional 
increase in Mach number that could have been observed in 
the absence of any correlations. 
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